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Abstract 
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the framework of the Lifshitz formula combined with the Drude dielectric function. An earlier result 
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I. INTRODUCTION 

During the last few years the Casimir effect has attracted a lot of experimental and 
theoretical attention as a nontrivial macroscopic evidence fo_r_ tlie ^existence of zer{>point 

and reviews 



)• 



oscillations of electromagnetic field (see, e.g., monographs 
Except for calculations of the Casimir force between perfectly shaped bodies made of ideal 
metal at zero temperature, the topical studies were conducted taking into account the effects 
of surface roughness, finite conductivity of the boundary metal and nonzero temperature . 



The correct theoretical description of the thermal Casimir force between real metals 



13 



assumed great importance due to recent precision experiments [tI-IsI ^ 10. 11. 1^ 
Q], followed by the prospective applications of the Casimir effect in nanotechnology 

3,0, 



has 
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and also its use as a test for predictions of fundamental physical theories 
13] • It was found unexpectedly that the calculations of the thermal Casimir force between 
two parallel plates made of real metal based on the Lifshitz formula ;24] supplemented 
by some model dielectric function ran into serious difficulties. The key question of the 
controversy is whether the transverse electric zero mode contributes to the Casimir effect 



in the case of real metals. In Refs. 



25 



26[ a positive answer to this question was obtained 



by the substitution of the plasma dielectric function into the Lifshitz formula. In the limit 



of infinitely high conductivity the results of Ref. j25|, |26|] are smoothly transformed into the 



familiar results for ideal metals j^j. In Ref. 



the Drude dielectric function was used to 



calculate the contributions of both longitudinal and transverse modes to the Casimir force. 
It was found that the transverse electric zero mode does not contribute 0|. The results of 
Ref. (2^ do not smoothly join with those for ideal metals. In the high-temperature limit, 
where the contribution of the zero modes is dominant, the Casimir force between plates 
made of Drude metals proves to be equal to one half of the force between the ideal metals. 
Today, there is an extensive literature on these two and other app roaches to the calculation 



of the thermal Casimir force between real metals l2f 



Q,H 3,QQ- 



HQ, 22,HQ,HHQ,HQ,Q, 



From this discussion it has been shown that the substitution of the Drude dielectric 
function into the Lifshitz formula results in negative values of the Casimir entropy within 
wide temperature interval and leads to the violation of the third law of thermodynamics 
(the Nernst heat theorem). On the contrary, Refs. have thrown doubt on the 



2 



computations of Ref. [39| and presented the numerical computations in favor of the statement 
that the Casimir entropy for the Drude metals is zero at zero temperature, i.e. the third 
law of thermodynamics is not violated. 

In the present paper we give a detailed and rigorous analytical derivation of the low- 
temperature behavior of the Casimir entropy in the framework of the Lifshitz formula com- 



bined with 
in Refs. 31 



;he Drude dielectric function (Sec. II). This derivation removes the doubts raised 



32 



and validates the thermodynamic inconsistency of the Drude dielectric 
function with the Lifshitz formula at nonzero temperature. The reason why the opposite 
conclusion is obtained in Refs. 0, 0, 0| is explained in Sec. III. This section contains also 



the resolution of the above thermodynamic puzzle 



jy presenting several computations in the 
Is? ! as opposed to the use of the Drude 



framework of the surface impedance approach |36 
model. The exact boundary conditions in terms of impedances, depending on polarization 
and angle of incidence, are compared with the Leontovich boundary condition. The latter 
is shown to be applicable to the case of fluctuating fields being in agreement with the third 
law of thermodynamics. The temperature dependences of the Casimir force and entropy for 
real metals in both approaches are compared with those for dielectrics. Sec. IV contains 
conclusions and discussion of recent experimental results. 



II. CASIMIR FREE ENERGY AND ENTROPY IN THE LIFSHITZ THEORY 
COMBINED WITH THE DRUDE MODEL 

The Casimir free energy for the configuration of two parallel plates at a separation a and 



temperature T is given by the Lifshitz formula 



2J|, which can be represented in the form 



J^{a,T) = ^Yl / ydy{ln[l-rl{Q,y)e-y]+\n[l-rl{Ci,y)e-y]}. (1) 
^'^ /=o "^fi 

Here prime adds a multiple 1/2 near the term with I = 0, is the Boltzmann constant, 
the dimensionless Matsubara frequencies are (i = C,i/uJc, where = 2'7ikBTl/h, Uc = c/(2a), 
and the reflection coefficients for the two different polarizations are expressed in terms of 
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the dielectric permittivity e{uj) as 

rliCi,y) = 



y- Vmi) + y' 



(2) 



We consider the Drude metals described by the dielectric permittivity of the Drude model. 
At the imaginary Matsubara frequencies it is given by 



1 + 



(3) 



0[0+7(T)]' 

where Up and 7(7) are the dimensionless plasma frequency and relaxation parameter defined 
by Up = Up/ujc, liT) = 'y(T)/uc. In the absence of relaxation 7(T) = and e^^'^ coincides 
with the dielectric permittivity of the free electron plasma model e^^\ Substituting Eq. 
into Eqs. (Q), (j21) one obtains the Casimir free energy J-'^^\a, T) and the reflection coefficients 
rj|'^^((^;, y) in the framework of the Drude model. If, from the very beginning, 7(T) = 0, then 
the Casimir free energy in the framework of the plasma model J-'^^\a^T) and reflection 



coefficients r^l']_{Ci,y) are obtained. It is notable that r^^'{0,y) = 0, whereas 



(o,z/) 



Spy 



iCuy) 



y - + y^ 



7^0, 



(4) 



y^ + y 

i.e. there is no smooth transition from T^^^ to JF^?') when 7(T) 0. This nonanalyticity 
is determined exclusively by the zero-frequency contribution of the transverse electric mode 
to the free energy T^^\a,T). 

For the calculation of the Drude free energy J-'^^\a,T) it is useful to represent it as 
the plasma free energy J-'^p'> (a, T) plus some additional terms. Taking into account that 



(0, y) = and (0, y) 



SDY 



{0,y) = 1, this identical representation is as follows 



^(^)(a,T) =^(p)(a,T) 



ydy In 



ApY 



io,y)e- 



(5) 



7 rr °° 



Svra^ 



1=1 



ydy \ In 



l-rf^\c,uy)e-y 



-In 



SpY 



{Chy)e~ 



+ ln \^l-rf ^\Q,y)e-y 



In 



ApY 



iChy)e 



An important point is that the zero-frequency contributions are contained only in the first 
two terms of the right-hand side of Eq. (jSJ, whilst the summation starts from the term 



with / = 1 leading to nonzero lower integration limits Q in the integrals at any nonzero 
temperature. 

It is instructive to find the asymptotic representation for the free energy applicable at 
T —>■ 0. First we notice that among the three parameters, contained in Eq. (jH)), i.e Up, (i (with 
/ > 1) and 7(T), the latter is the smallest one. In fact at T = 300 K for good metals 7 ~ 
1013 _ iQi^rad/s (for Au 7 = 5.32 x lO^^rad/s) whereas ^1 = 2nkBT/h = 2.46 x lO^^rad/s, 
and C,i = Z^i, i.e. in all cases 7 < When T decreases from room temperature up to 
approximately Td/4, where is the Debye temperature (Tp = 165 K for Au jlll), 7(T) ~ 
T, i.e. 7(T) decreases following the same law as ^i, preserving the inequality 7(r) < At 
T < Td/4 the relaxation parameter decreases even more quickly than C,i with decreasing T 



(as ~ according to the Bloch-Griineisen law due to electron-phonon collisions 



~ at liquid helium temperatures due to electron-electron scattering j4l|). As a result, at 
T = 30K 7(T)/^i(T) ^ 4.9 x IQ-"^, at T = 10 K 7(r)/^i(r) ^ 1.8 x lO'^, and this relation 
decreases further with T — 0, i.e. at low temperatures the condition 7(T) ^ ^i(T) is largely 
satisfied. 

The largest parameter of the above three is Up (for Au Up = 1.37 x lO^^rad/s). For 
example, for Au at T = 300 K, 70 K, and 10 K we have 7(T)/cUp = ^{T)/ujp = 3.88 x 10~^ 
6.71 X 10-^ 1.06 X 10-^ respectively. 

The refiection coefficients (j21) have continuous derivatives with respect to the relation 
7(T)/,^i(T) (/ > 1) at the point 7(T)/^i(T) = 0. Under the above proved condition 7(T) <^ 
^i(T), which is satisfied at all sufficiently low temperatures, the relation 'y(T)/^i(T) -C 1, 
and we can expand in Taylor series around a point 'y{T)/C^i{T) = keeping only the first 
order terms 



42 1 and as 



f\Cuy)-^^R\\{Cuy). 
Yi(i,y)~^RdCi,y), 



where 



R\\iCi,y) 



R±{Ci,y) 



CiiTY 

2Cfay[l + a^2y'-Cf)]\rP{Q,y)\ 



(6) 



y + aQ [oiy + a/1 + Oi^y"^ 



2 ' 



2«y|rf (0: 



a/1 + a^y"^ 



{ay+./l + 



a^y"^ 



2 ' 



(7) 



and a = l/cDp = Ap/(47ra), \p is the plasma wavelength (note that we keep the argument T 
when (i participates in the expansion parameter but omit it in the other cases). 
The same expansions for the logarithms which appear in Eq. (Q) are 

7(r) R\\{Ci,y)e-y 



In 



In 



1 - {Chy)e 



l-rf^ {Ci,y)e~ 



In 



In 



{Ci,y)e~ 



+ 



+ 



CK^)l-ri^)'(0,y)e-^' 

7(T) R±{Ci,y)e-y 



(8) 



Substituting Eq. (jSI) in Eq. © we obtain the following expression for the Casimir free 



energy between parallel plates made of Drude metal 



ydy In 1 - rf\o,y)e~y + J^^^\a,T), (9) 



where the contribution depending on the relaxation parameter is given by 



ydy 



R\\iCi,y) 



+ 



ey-rP\Ci,y) ey-rf\ci,y) 



(10) 



Notice that in this expression the small parameter 'y{T)/ujp = '~f{T)/ujp, which does not 
depend on an index of summation, is put in evidence. 

The low-temperature asymptotic limit of the Casimir free energy J-'^^\a, T) in the frame- 
work of the plasma model was investigated with details in Ref. j^Sf (the coinciding numerical 



results follow also from the computations of Ref. |25|). In terms of the two small parameters 
T/Teff [where ksTeff = hc/{2a) = fkVc] and 6o/a = 2a [where 6o = Ap/(27r) is the skin 
,aye. thickne. .n the frequency .epon of the i,A„ed optics) the pla.na .nodel fee energy 
is given by |39|| 
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^cC(3) 



1 + 2 



5o 



T 



eff 



TT 



45C(3) 



1 + 4- 



.5n 



T 



eff 



, (11) 



where E^^^a) is the Casimir energy at zero temperature calculated by using the plasma 
model dielectric function. The perturbation expansion (fTTj) is applicable at separations 
Ap < a < 3 jjia at all temperatures T < 300 K |2fil |. 

The second term in the right-hand side of Eq. Q is linear in the temperature. It is 
an easy matter to calculate the coefficient near T perturbatively. For this purpose we use 



6 



Eq.(0]) and expand the logarithm under the integral in powers of 60/ a. Then all integrals 
are taken explicitly, resulting in 

2 



ydyhi 



(0,l/)e- 



W(3) 



.5n 



41:1 + 12 



5n 



-32 ('^ 



C(5) 
16C(3), 



+ 80 ('^ 



C(5) 
4C(3) 



where is the Riemann zeta function. 

Now we consider the low-temperature behavior of the last term in the right-hand side of 
Eq. 0, J-'^'^\a,T), which depends on the relaxation parameter. An important point is that 
the relaxation parameter is not an independent one but is the function of the temperature, 
7 = 7(T). Because of this, it is not improbable that the temperature dependent terms 
resulting from the integration and summation in Eq. (jlUj) will cancel the second term, linear 
in the temperature, in the Casimir free energy 0. Below we demonstrate that this is not 
the case. 

As pointed out above, 7(T) —>■ when T ^ no slower than ~ T^. If desired that the 
quantity J-'^^^ (a, T) from Eq. (fTUj) be linear in T, the sum in Eq. ()10j) should tend to infinity 
as when T — > 0. Let us find what is the actual asymptotic behavior of the quantity 
J^(T)(a,T) when T 0. For this purpose we expand R\\^_i{(i,y) from Eq. ((7j) up to the 
first order in the small parameter a (recall that J-'^^^ is already proportional to the smallest 
parameter of our problem 'y/ujp) 

2Qa 



R\\{Cuy) 



y 



Substituting Eq. (fT^ into Eq. (fTUj) one obtains 



P 1=1 



EC 



R±{Ci,y) = 2?/Q!. 



dy 



(13) 



y^dy 
ey -1 



(14) 



Each of the two sums in Eq. (fT^ can be simply found when T — > 0. The asymptotic of the 
first sum is as follows 



1=1 



dy 



2tiT 



ey -I T, 



00 ^ 



k=l 



+ 



1 



27rT 



C(3)+C(2) (15) 



(here we have neglected all terms in the expansion of exp{2'n'kT/Teff) in the denominators 



7 



starting from the third ones). For the second sum in Eg . (1141) one obtains 



1=1 



y^dy 
ey -1 



eff 



27tT 



^ P ^ ^ Teff ^ k' 

k=l k=l 



^2 ^2^kT/T,fs _ I 



(16) 



+ 



Teffm j^ 27rT 



3Te//C(3) 



vrT 



TkT/T,ff _ 

Substituting Eqs. (fTH|l . (fTH|) into Eq. (fT^ we arrive at 
:(7).. ^^ ^ 7(T) A:^r,^/C(3) 



eff 



LUr, 



, 27rT 
- In + 2 



37r 



27rT 



C(2) T 



+ 2C(2). 



(17) 



As is seen from Eq. (fT7j) . the leading term in J-'^"'\a,T) behaves as — 7(T) ln(T/Tejj) and 
goes to zero when T — >■ because 7(r) ~ at hehum and lower temperatures. 

Now we are in a position to find the low temperature behavior of the Casimir entropy 

9^(^)(a,T) 



S^''\a,T) 



dT 



(18) 



calculated by using the Drude dielectric function, where the Casimir free energy JF^^) (a, T) 
is given by Eq. 0. According to Eq. 0, there are three contributions into the Casimir 
entropy in the framework of the Drude model at low temperatures. The first one is given 
by the Casimir entropy calculated by means of the free electron plasma dielectric function. 
It is obtained from Eq. (Illj) 

aj^(f)(a,T) 3kBC{3) [ T 



S^P\a,T) 



X 



dT 
47r3 T 



87ra2 



eff 



(19) 



.5n 



87r=^ 



T 



135C(3)T,;yJ 



39| (it is also in agreement with the 



25, 



This result coincides with the one obtained in Ref. 

numerical computations of the thermal corrections to the Casimir force in Refs. 
Evidently the plasma model Casimir entropy is positive and S^P\a,T) when T — 0, 
i.e. it is in agreement with the Nernst heat theorem. 

The second contribution to the Drude model Casimir entropy is obtained from the second 
term in the right-hand side of Eq. © taking into account Eq. (jl2j) and is given by 



5r(a,T) = 5r(a,0) 



k 



B 



MM 



1-4- 



12 



ydyln 

2 



l-r'f^'{0,y)e-y 



(20) 
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[see Eq. ffT^ for higher perturbation orders in So /a]. This contribution is negative and does 
not depend on the temperature. 

The asymptotic behavior of the last contribution to the Casimir entropy in the Drude 
model, as given by Eq. ©, at T — * 0, is obtained from Eq. (jl7|) with account 7(T) = 7oT^ 



, 2-kT C(2) T 

-21n-— + 3 + 97r>^ ' 



Teff C(3) T, 



eff 



21) 



dT 47r2a2 ujp T 

From Eq. (PT|) we notice that S^"'\a,T) when T ^ 0. 

As a result, the value of the Casimir entropy at zero temperature calculated with the 
help of the Drude dielectric function is found from Eqs. (fT^ - (PT|) : 

(a, 0) = lirn^ [S^P^ (a, T) + 5^""^ (a, T) + 5^^^ (a, T)] = 5^""^ (a, 0) < 0, (22) 

where the quantity in the right-hand side is negative and is given by Eq. ((201) • This quantity 
depends on the parameters of the system, such as the separation between the plates and 
plasma frequency, violating the third law of thermodynamics (the Nernst heat the- 

orem). Therefore we may conclude that the Drude dielectric function is thermodynamically 
inconsistent with the Lifshitz formula and unusable to calculate the thermal Casimir force 
between real metals. Note also that according to Eqs. ()20 |) -()22 |) the Casimir entropy be- 
tween two parallel plates made of Drude metals is negative within a wide separation range. 
The entropy cannot be made positive by introducing some composite system containing a 
subsystem with negative entropy between the two infinite plates without changing the value 
of the Casimir force. 

In the next section we discuss the possibilities to avoid the above thermodynamical puzzle 
and formulate the approach in which the Casimir entropy is positive, becoming zero at zero 
temperature. 



III. TOWARDS A THERMODYNAMICALLY CONSISTENT THEORY OF THE 
THERMAL CASIMIR FORCE BETWEEN REAL METALS 

The above asymptotic limit for the Casimir entropy was derived under the important 
condition that 7(T) <^ ^i{T), i.e. at any T the magnitude of the relaxation parameter must 
be much less than the Matsubara frequencies. If this condition does not hold, the obtained 
conclusion concerning the thermodynamical inconsistency of the Drude model combined 

9 



with the Lifshitz formula is open to question. This explains why in Refs. 



31 



3 it 



32, 3J it was 



concluded that the Lifshitz formula combined with the Drude model respects the Nernst heat 



theorem. In Refs. 



31 



32, 



33| all computations were performed under the condition that the 



relaxation parameter is constant (the values 7 = 35.6 meV or 0.01 meV, as at T = 300 K 
and T = lOK, respectively, were extended to all temperatures). Then even for the smaller 
value, at temperatures T < 0.018 K the inequality 7 > ^i(T) is fulfilled, i.e. the condition 
7(T) ^ C,iiT) is violated, and our proof in Sec. II is not apphcable. 

The question arises whether there are physical prerequisites for violating the condition 
7(T) <^ C,i(T). According to Sec. II, the nonelastic processes of electron-phonon collisions 
and also the elastic electron-electron scattering respect the inequality 7(r) ^ ^i{T)- One 
may hope, however, that some fine properties of real metal bodies could result in the violation 
of this inequality at some sufficiently low temperatures. It has been proposed j3] that this 
role can be played by impurities and defects which lead to a nonzero residual value of the 
static resistivity (and, thus, a relaxation parameter; recall that resistivity is proportional 
to the relaxation parameter as the temperature goes to zero (see also Appendix D in 
Ref. 0)- 

Here we adduce the argument that impurities cannot remedy the situation with the 
violation of the Nernst heat theorem. In fact, the resistivity ratio of a sample can be 
defined as the ratio of its resistivity at room temperature to its residual resistivity. For pure 
samples the resistivity ratio may be as high as 10^ [41]. As an example, let us consider 
Au with 7(T = 300 K) = 5.32 x lO^^rad/s. In this case for the residual value of the 
relaxation parameter one obtains 7r.es = 5.32 x lO^rad/s. The asymptotic expressions of 
Sec. II are applicable under the condition 7 <^ ^1 = 2TTkBT/h. Thus, with allowance 
made for impurities, these asymptotics are applicable at temperatures T ^ h'-jres/ {'^'^ks) = 
6.5 X 10~^K. What this means is that the Casimir entropy at temperature T = 5 x 10~^K 
has a nonzero negative value given by Eq. ()20|1 . Physically this is equivalent to the violation 
of the Nernst heat theorem. Only at smaller temperatures of about 10^^ K does Eq. 
break down and the Casimir entropy rapidly takes a sharp upward turn to zero. We would 
like to point out also that the usual theory of the electron-phonon interaction, describing 
electrons interacting with elementary excitations of a perfect lattice with no impurities, must 
satisfy and does satisfy all the requirements of thermodynamics. In fact, the mean free path 
of the electrons between the collisions with impurities at zero temperature is many orders 
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of magnitude greater than the penetration depth of the electromagnetic oscillations at the 
characteristic frequency into the metal. To make sure that this is the case, recall that the 
relaxation time for Au (the inverse of the relaxation parameter) at T = 300 K is equal to 
t(T = 300 K) = 1.88 X lO^^^s. Using the above resistivity ratio, one obtains the value of 
the relaxation time at zero temperature r(T K) = 1.88 x 10~^ s. Finally, using the value 
of Fermi velocity for Au {vp = 1.79 x lO^m/s) the mean free path of the electron is equal to 
1{T ^ OK) = vft{T ^ ok) = 3.36cm. This should be compared with the thickness of the 
skin layer in the region of the infrared optics, equal to approximately 22 nm. That is why 
the attempt to remedy the violation of the Nernst heat theorem at the expense of impurities 
is meaningless. 

Recently the resolution of these complicated problems was obtained |37j using another 
approach to the description of real metals based on the concept of the Leontovich impedance 
boundary conditions. This approach offers a fundamental understanding of the reason why 
the Drude model is not compatible with the theory of the thermal Casimir force between 
real metals. 

The main concept of the Lifshitz theory is the fluctuating electromagnetic field consid- 
ered on the background of dielectric permittivity depending only on frequency. This concept 
works good in the case of dielectrics but is not adequate for real metals. In fact, in the fre- 
quency region of the anomalous skin effect the spatial non-uniformity of the field makes 
impossible a description of a metal in terms of ^(co') Then the electromagnetic fluc- 

tuations also cannot be considered on this background. Moreover, in the frequency region 
of the normal skin effect the electric field E initiates a real current of conduction electrons 
j= ctqE, where ctq is the DC conductivity of a metal. These j and E should be considered 
as real ones |32]- In contrast with E, the fluctuating field cannot heat a metal as E does 
due to collisions of conduction electrons with phonons. 

At present a complete theory of field quantization inside metals which, among other 
things, should take into account the effects of spatial non-uniformity, is not available. Be- 
cause of this, the concept of the electromagnetic fluctuations inside a metal remains unclear. 
In the absence of a complete theory we should not take into consideration the metal inte- 
rior, but rather take into account the realistic material properties by means of the surface 
impedance function. 

It is well known that for a plane wave of a single frequency inside a medium with dielectric 



11 



permittivity e the following equations are valid 

ujH = ckx E, coeE = -ck x H, (23) 



where fc is a complex wave vector. Then from the first equality of Eq. ()23|) for the field with 
transverse polarization {E is perpendicular to the x2;-plane which is the plane of incidence) 
inside a metal near its boundary plane it follows 

E, = [H, X n] , Z^= (24) 

Here the index t refers to the component of the field parallel to the boundary plane and the 
unit vector n is perpendicular to it and directed inside a metal. In the same way, from the 
second equality of Eq. (j23|) for the field with longitudinal polarization one obtains 



The quantities Z^^\\ are called impedances 
identically represented as 



3- 



n X Et] = Z\\Ht, Z\\ = —Ju^e - c^kl. (25) 



By the use of the Snell's law they can be 



1 1 / sin^?9o . ^ 

^ll = ^Vl 26 



where -t^o is the angle of incidence of the electromagnetic wave from vacuum on the boundary 
plane of the metal. 

In metals for all frequencies which are at least several times less than the plasma frequency 
we have that \e\ ^ 1. For this reason, the term sin^-i^o/^ in Eqs- (ESI) can be neglected in 
comparison with unity. This leads to the fact that inside a metal all waves are spreaded 
perpendicular to the surface, i.e. the refraction angle is equal to zero independently of the 
angle of incidence 

Et = Z [Ht X n] , [n X E^] = ZH^ (27) 
with Z = 1/ ^Je can be used as boundary conditions in order to determine the field outside 



the metal. The quantity Z = Z{uj) is called the impedance of a metal or the "intrinsic" 
impedance 4^1 ■ It depends only on u and does not depend on the polarization or the angle 
of incidence. We emphasize that for real photons the difference between the Leontovich 
impedance, as is in Eq. and impedances in Eq. is negligibly small. What is more, 

12 



when cij — > 0, the dielectric permittivity goes to infinity and the Leontovich impedance 
coincides precisely with the impedances (j^ . By postulating the boundary condition ^T7\ 
in the theory of the Casimir effect, we admit, in fact, that the virtual photons have the same 
reflection properties on the metal boundary as real ones do. It is signiflcant that the surface 
impedance and the boundary conditions (|27p still hold, even in the frequency domain of the 
anomalous skin effect when, due to the spatial non-uniformity of the fleld, the description 
in terms of e becomes impossible. Notice that in recent Ref. ||4^ it has been suggested 
to use the "exact" boundary conditions (j211), with the impedances depending on the 
wave vector (he. on the angle of incidence) instead of the Leontovich conditions (j^7j) used 



in Refs. 



• 3 



37|. This, however, leads us back to all the above problems with the thermal 



Casimir force, because in the absence of a mass-shell equation the representation for the 
impedances (j^ . (^3)) becomes not equivalent to the representation (pUjl . In our case we 
have used the Leontovich boundary condition (j^Tj); that is to say, the representation (PU)) 
was generalized for the case of virtual photons. If one generalizes Eqs. fl24j) . ()25|1 for the case 
of virtual photons, this will lead to zero value for the transverse reflection coefficient at zero 
frequency and, therefore, to a contradiction with thermodynamics. 

By the use of the surface impedance instead of the Drude model 0, the Lif shitz formula 
(HJ is preserved, but the coefficients r^, r'j_ given by Eq. © should be replaced by 



y - zmci 

[y + Z{zOKi 



- zmy ' 

[Ci + Z{zCi)y\ 



(28) 



Substituting in Eq. (j^Hj) the impedance function of the normal skin effect or the anomalous 



skin effect, one finds r?:{0,y) = r'j_{0,y) 



follows 



1 0, 



37|. In the region of the infrared optics it 



that r2(0, y) = 1, ri(0, y) = {up - yf/{Cop + yf. 



It should be stressed that the expressions obtained for the refiection coefficients at zero 
frequency in the impedance approach are exact. They readily follow from the exact Eq. (pUj) 
as I sin^ j^n /g I < {"^ /uj'^)uj when a; — 0. This result is in contradiction with the statement 



of Ref. [33[ . The authors of Ref. [33[ start from representation ()24|) for Z± and consider the 
limit a; — >■ at fixed nonzero k± [i.e. they violate the mass-shell equation from which 
Eqs. (j^l), are equivalent to Eq. (jSH))]- As a result, the equality r^^O, k±) = obtained 
in Ref. 23] leads to the violation of the third law of thermodynamics (see Sec. II). In fact, 
both approaches, ours and that of Ref. js^, start from different postulates. Our postulate 
is that the reflection properties for the fluctuating fleld are the same as for real photons. If 
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this is true, the impedances ()26|) follow, which coincide with the Leontovich impedance at 
zero frequency and are approximately equal to it at all nonzero frequencies with a very high 
precision [all calculational results based on Eq. and the Leontovich inpedance using 
Z = \ j ^/e for metals are practically the same] . In Ref . jssj another postulate is assumed, 
which admits that the reflection properties for the fluctuating field are different from those 
for the usual electromagnetic field. Both postulates have the right of being assumed because 
it is impossible to study the reflection properties of the virtual photons experimentally. The 
second postulate, however, is shown to be inconsistent with thermodynamics. For this reason 
it must be rejected, and so we conclude that the fluctuating field has the same reflection 
properties on a metal boundary as the usual electromagnetic waves. 

Let us now present several computational results obtained using Eqs. (jT} and in the 
framework of the impedance approach in comparison with the results calculated from the 
Lifshitz formula (0), (j21) combined with the Drude model Q. In Fig. 1, the magnitude of 
the Casimir force per unit area F = —dT /da for gold plates is plotted, in the temperature 
interval 1 K< T < 1200 K at a separation distance a = 1 /xm. At this separation distance, the 
characteristic frequency ujc belongs to the region of the infrared optics where the impedance 
function is given by 

Z{iCi) = (29) 



The solid line represents the values calculated in the framework of the impedance approach, 
and the dashed line is obtained via the Lifshitz formula supplemented by the Drude model 
with a temperature dependent relaxation parameter [data for 7(T) are taken from Ref. 
It is clearly seen that the dashed line is not monotonous, demonstrating the existence of a 
wide temperature region where the force modulus decreases with an increase of the temper- 
ature (as in Figs. 2,3 of Ref. 3^). At the same time, the solid line obtained by using the 
surface impedance demonstrates the monotonous increase of the magnitude of the Casimir 
force with temperature in perfect agreement with what is expected from thermodynamics. 

It is instructive to know if the nonmonotonous force-temperature relation takes place 
only for Drude metals or if this may also happen for dielectrics. In Fig. 2, the magnitude 
of the Casimir force per unit area between dielectrics with e = const is shown for different 
temperatures at a separation of a = 1 fim. Both solid and dotted lines were obtained from 
the usual Lifshitz formula (Q), Q with e{i(i) = e = const. The solid line is for mica with 
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e = 7; the dotted line coincides with the hne in Fig. 5 of Ref. l33 with e = 100. The sohd hne 
shows a monotonous increase of the Casimir force with the temperature, as is expected from 
thermodynamics. On the other hand, the dotted hne is not monotonous as the dashed hne 
in Fig. 1. This is, however, artifactual because e can be assumed to be independent on the 
frequency and temperature only in the case of so-called non-polar dielectrics whose atoms 
or molecules do not have their own dipole moments. The electric susceptibility of non-polar 
dielectrics arises due to the electronic polarization of atoms and molecules. The values of 



49 



50|,|51|. Large values of e can exist 



e for non-polar dielectrics are of the order of one 
only for polar dielectrics where the partial orientation of permanent dipole moments occurs. 
But for polar dielectrics e depends strongly on the frequency and temperature. Specifically, 
their e quickly decreases with the increase of frequency. As a result, at optical and infrared 
frequencies, which are characteristic for the Casimir effect, the values of e are determined by 
the electronic polarizability j49j and cannot exceed several units. The dielectric permittivity 
of the polar dielectrics along the imaginary frequency axis can be modeled by 

^(^o = i+Er^' (30) 

n=l Uln 

where the sum describes the effect of possible Debye rotational relaxation frequencies (the 
absorption spectra of dielectrics are not influential at the characteristic frequencies of the 
Casimir effect uj^ when the separation between plates is of the order of 1 /im). In solids, the 
polarization due to orientation of the permanent dipole moments disappears at rather low 
frequencies ~ (10^ — 10^)rad/s At very high frequencies uo > lO^^rad/s, e decreases 
to unity 0]. Therefore, for the model calculation we may choose N = 2, Ci = 93, ui = 
lO'^rad/s, C2 = Q, 002 = lO^^rad/s. In this case e{0) = 100 and £ ~ 7 at the characteristic 
frequency Uc corresponding to a = 1 fim. 

Eq. (jHUj) was substituted into the Lifshitz formula ((H), In Fig. 2, the magnitude of 
the Casimir force per unit area between the polar dielectrics is represented by the dashed 
line. It is seen that the force-temperature relation, given by the dashed line, is monotonous 
as in the case of non-polar dielectrics with small e, as expected from thermodynamics. From 
Fig. 2 it is clear that the Casimir force between dielectrics is a monotonous function of the 
temperature when realistic input data for e are substituted into the Lifshitz formula. 

In Fig. 3, the Casimir entropy for gold is plotted as a function of the temperature at a 
separation distance between plates a = 1 fim. The solid line is drawn in accordance with the 
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impedance approach, i.e. by the use of Eqs. ((T)), (j2Hl) and ((221) • The dashed hne is obtained 
from the usual Lifshitz formula and Drude model (QJ-Q- Evidently, the solid line satisfies 
all conditions, i.e. positive values of the entropy at nonzero temperatures, and the validity 
of the Nernst heat theorem. By contrast, the dashed line presents the negative values of the 
entropy and the violation of the Nernst heat theorem. The analytical proof of the validity 
of the Nernst heat theorem in the impedance approach can be found in Ref. js?]. 

IV. CONCLUSIONS AND DISCUSSION 

As we have proved above, the substitution of the Drude dielectric function into the 
Lifshitz formula for the thermal Casimir force leads to the violation of the third law of 
thermodynamics (the Nernst heat theorem). A rigorous analytical evidence of this statement 
lies on the fact that at low temperatures the magnitude of the relaxation parameter is much 
less than the Matsubara frequencies, a property which is always true in the case of perfect 
crystal lattice. A special analysis of the role of defects or impurities leads to the conclusion 
that they are incapable to reconcile the calculations of the thermal Casimir force between 
Drude metals with thermodynamics. 

It has been known that the Lifshitz formula combined with the Drude dielectric function 
predicts a linear (in temperature) large thermal correction to the Casimir force at short 
separations connected with the second term in the free energy 

recent precision measurement of the Casimir pressure between Au and Cu plates at room 
temperature jisl. llG^ this correction (which we call "alternative" ^^) comprises 4.89 mPa 
and 1.23 mPa at separations a = 300 nm and 500 nm, respectively, i.e. 3.6% and 7.24%, 
respectively, of the total Casimir pressure. For reference, the traditional thermal corrections 
in the same experimental configuration (computed using the plasma dielectric function [2^ 
13] or the surface impedance 0, Q) at separations a = 300 nm and 500 nm are equal to 
-0.00863 mPa and -0.00441 mPa, respectively, i.e. only 0.006% and 0.03%, respectively, of 
the total Casimir pressure (note that the traditional thermal corrections are of the same 
sign as the Casimir force, i.e. negative). The experimental results of Refs. jisl 16] and the 
extent of their agreement with theory rule out the linear thermal correction to the Casimir 
force as predicted by the Drude dielectric function substituted into the Lifshitz formula (see 
Ref. for details). Thus, the combination of the Drude model with the Lifshitz formula 



16 



is not only thermodynamically inconsistent, but is also in contradiction with experiment. 

As we have demonstrated through several computations, the description of real metals 
based on the surface impedance (Leontovich boundary conditions) is thermodynamically 
consistent. Unlike the dielectric permittivity, the surface impedance is well defined at all 
frequencies, even in the domain of the anomalous skin effect. By the use of the Leontovich 
impedance boundary conditions, instead of the Drude model, the Lifshitz formula is pre- 
served, but the reflection coefficients are expressed in terms of the impedance function. An 
important point is that different representations for the impedances, which are equivalent 
for real photons, become nonequivalent in application to the fluctuation fields. As we have 
shown above, in the case of fluctuating fields the representation (|26|) for the impedance 
should be used leading to the Leontovich boundary conditions, whereas the representations 
(j211), when applied to virtual photons, lead to the contradictions with thermodynamics. 

In the framework of the impedance approach the value of the zero-frequency term of the 
Lifshitz formula is prescribed by the form of the impedance and quite satisfactory physical 
results are obtained. In particular, the Casimir energy and force turn out to be monotonous 
functions of the temperature, in agreement with what would be expected from thermo- 
dynamics. In the region of the infrared optics the results obtained by using the surface 
impedance coincide with those found earlier through the use of the plasma dielectric func- 



tion 
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The Casimir entropy in the impedance approach is always positive and vanishes 
at zero temperature, in accordance with the Nernst heat theorem. 

To conclude, we have proved that the Drude dielectric function is not appropriate to 
describe the thermal Casimir effect in the case of real metals, leading to contradictions with 
thermodynamics. On the other hand, the Lifshitz formula with the coefficients expressed in 
terms of the surface impedance is suitable to calculate all the quantities of physical interest. 
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Figures 
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FIG. 1: Magnitude of the Casimir force per unit area for gold versus temperature when a = 1 /im. 
The sohd line is obtained in the framework of the impedance approach. The dashed line is obtained 
by using the Drude dielectric function. 
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FIG. 2: Magnitude of the Casimir force per unit area for dielectrics versus temperature when 
a = 1 iJ,m, calculated by the Lifshitz formula. The solid line is for mica, with e = 7. The dotted 
line is obtained for non-existent non-polar dielectric with e = 100. The dashed line is for polar 
dielectric with £(0) = 100. 
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FIG. 3: Casimir entropy per unit area for two gold semispaces versus temperature, when a = 1 /xm. 
The sohd hne is obtained in the framework of the impedance approach. The dashed hne is obtained 
by using the Drude dielectric function. 
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